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ABSTRACT 

We present the results of polarimetric (R band) and multicolor photomet- 
ric (BVRIJHK) observations of the blazar AO 0235+16 during an outburst in 
2006 December. The data reveal a short timescale of variability (several hours), 
which increases from optical to near-IR wavelengths; even shorter variations are 
detected in polarization. The flux density correlates with the degree of polariza- 
tion, and at maximum degree of polarization the electric vector tends to align 
with the parsec-scale jet direction. We find that a variable component with a 
steady power-law spectral energy distribution and very high optical polarization 
(30-50%) is responsible for the variability. We interpret these properties of the 
blazar within a model of a transverse shock propagating down the jet. In this 
case a small change in the viewing angle of the jet, by < 1°, and a decrease in 
the shocked plasma compression by a factor of ~1.5 are sufficient to account for 
the variability. 

Subject headings: galaxies: active — galaxies: BL Lacertae objects: individual 
(AO 0235+164) — galaxies: jet — polarization 



1. Introduction 

Blazars are active galactic nuclei (AGNs) that possess extreme properties such as vio- 
lent variability, high polarization, relativistic jets, and bright 7-ray emission. It is commonly 
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thought that the radio through optical emission in blazars is synchrotron radiation that 
originates in the jets. Polarimetric observations provide crucial data for investigating the 
magnetic field structure that forms the framework of such radiation. However, the polar- 
ization behavior is usually too complex for straightforward interpretation. For example, 
long-term observations of BL Lac reveal a change of degree of polarization from 0.5% to 
>20% with all possible electric vector position angles (EVPAs), although the distribution 
of EVPAs indicates the existe nce of a preferential direction aligned with the parsec-scale 
jet (IHagen-Thorn et al. Ill994l ). Difficulties are even more severe for physical models aimed 
to explain both the brightness and polarization behavior. A connection between polariza- 
tion and brightness is rarely seen at optical wavelengths (jSmith I Il996l ). However, there 
are strong indicatio ns that the optical po l arized emission is linked to the mm- wave VLBI 



core of blazars jets (ILister fc Smith 1 12000| ; iGabuzda et al. I l2006t iD'Arcangelo et al. 
Jorstad et aDl2007t ). 
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At radio wavelengths some progress in interpetation of both polarization and brightness 
variability has been achieved by using a model in which a shock propagates down the jet (e.g., 
Hughes. Aller. fc Aller Ill989l ) . At optica l wavelengths, where polarization and brightness can 



vary on timescale of a few hours (e.g., lHagen- Thorn Ill972l ; iMiller fc Noble Ill996l ; ICarini 



20061 ). a proper test of models requires very intensive monitoring of both polarization and 
brightness. Such simultaneous photometric and polarimetric measurements are difficult to 
obtain, but they provide a vital tool to distinguish among different models responsible for 
the variability. 

We have observed the BL Lac object AO 0235+164 (z=0.94) during the relatively 
short outburst in 2006 December as part of the polarimetric and photometric monitor- 
ing program of blazars that has been carried out at Saint Petersburg State University for 



many year s. AO 0235+ 



ity (e.g., IRaiteri et al. 
(26 ± 7)c jPiner et al. 



64 s hows extreme blazar properties, such as intraday variabil- 
2005h and very high apparent speeds in the radio jet, /3 app = 
20061. that have drawn the attention of many investigators to 



this o bject (IRaiteri et al. Il200ll . I2005L 120061 : ljunkkarinen et al. 1 12004 ; lHagen- Thorn et al 



20071 ). In addition, AO 0235+164 is one of the most highly polarized BL Lac objects known. 
During the outburst in 1983 Jan uary the opti cal polarization of the object varied from 12% 
to 24% over a 24 hour period (jSmith I Il996l ). Our polarimetric monitoring program has 
caught the source in a similar state but with an even higher degree of polarization, > 30%. 
The polarimetric observations are enhanced by multicolor photometric measurements. Here 
we present the observations and show that the correlations that we have uncovered provide 
strong evidence in support of the shock model. 
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2. Observational data 

Observational data at optical wavelengths (BVRI bands) were obtained at the 70- 
cm reflector of the Crimean Astrophysical Observatory with the photometer-polarimeter 
of the Astronomical Institute at St. Petersburg State University based on an ST-7 CCD. 
M ethods of observation and d ata reduction of the photometric mode are described in detail 
Hagen-Thorn et al. I (bood ). We used the comparison stars recommended by the Whole 



in 



Earth Blazar Telescope (WEBT) program 0. The photometric errors do not exceed 0.03 r 
B band and 0.02 m in other bands when the object is brighter than 17 m . 

Observations at near-IR wavelengths (JHK bands) were carried out at the 1.1-m tele- 
scope of the Main (Pulkovo) Astronomical Observatory of the Russian Academy of Sciences 
located at Campo Imperatore (Italy) and equipped with the SWIRCAM camera. Meth- 
ods of reduction are the same as for the optical observations. Errors in JHK do not ex- 
ce ed 0.02 m when the object is brighter than 16 m in J band. In all cases the calibration 



of iMead et al. I (Il990l ) was used to transform magnitudes to flux densities (subsequently 



referred to as "fluxes"). 

Polarimetric observations were performed in R band using two Savart plates rotated 
by 45° relative to each another. By swapping the plates, the observer can obtain either 
the relative Stokes q or u parameter from the two split images of each source in the field. 
Instrumental polarization was found via stars located near the object under the assumption 
that their radiation is unpolarized. The Galactic latitude of AO 0235+164 is —39° so that 
interstellar polarization (ISP) in this direction is less than 0.7%. To correct for the ISP, the 
mean relative Stokes parameters of nearby stars were subtracted from the relative Stokes 
parameters of the object. This accounts for the instrumental polarization as well. The errors 
in degree of polarization are less than 2% for R < 16.5 m , while the EVPA is determined with 
a precision of several degrees. 

The results of photometric observations in the optical and IR bands are given in Tables 
Q] and [21 respectively. The light curves are shown in Figure [T] for the convenience of visual 
inspection. One can see that during the outburst the brightness increases by more than 
3 m . Table |3] lists the results of the polarimetric observations. In Figure [2] the upper panels 
illustrate the photometric R and polarimetric observations for three consecutive Julian dates 
during the most dramatic period of the outburst. Intraday variations (IDV) are evident 
within a given Julian date. The bottom panels present absolute Stokes parameters during 
the same period. The existence of systematic variations during each night is a prominent 



1 http://www.to.astro/blazars/webt/ 
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feature of all the curves. 



3. Global Properties of Variability 
3.1. Timescale of Variability 

We have determi ned the timescale of variability o f the flux in each band using the 



definition proposed by lBurbidge. Jones, fc O'Delll (119741 ): r = rft/ln(F 1 /F 2 ), where dt is the 



time interval between flux measurements F\ and F 2 , with Fi > F 2 . We have calculated 
all possible timescales 7y for any pair of observations for which |Fj — Fj\ > + &Fj at 
frequency v. The minimum timescale at frequency v is determined as r„ = minjrij^}, where 
i—l,...N — 1, j—i+l,...N, and N is the number of observations. Table H] lists the results of 
these calculations as well as the timescale of the polarized flux variability. The uncertainties 
of the values of r v are determined by errors of the flux measurements. However, the minimum 
timescale of variability that we derive is mainly constrained by the sampling, for example, 
the B, V, and / light curves have only two measurements during JD 2454077 (Tabled]), the 
date at which we find the minimum timescales in R band and at IR wavelengths. In Table 
H] we enclose the timescales derived from these light curves in brackets to indicate that they 
can not be compared with the timescales derived from the better sampled R, J, H, and K 
light curves. 

The values of r v indicate that during the event the timescale of variability was only a 
few hours and the shortest timescale, ~0.3 hr, occured in the polarized flux. Comparison of 
the timescales of the total and polarized flux variability in R band implies that the polarized 
emission originates in a subsection (~l/5) of the optical emission region. Moreover, the light 
curves with similar sampling during the event (R, J, H, and K data) reveal an increase of 
timescale with wavelength. 



3.2. Multicolor Behavior 

We have calculated the discrete cross-correlation function, DCF, between the R band 



l ight c urve and those at other wavelengths using the approach suggested by lEdelson fc Krolik 



(119881 ) for unequally spaced data. The functions were computed for lags between light curves 
from —6 to +6 hr with a step of 0.5 hr and with the constraint that at least 5 pairs of measure- 
ments are needed to define the coefficient of correlation at a given delay. Figure [3] presents 
the results of the calculations. The absence of points at some lags indicates that the data 
are not capable of determining the DCF for such delays owing to the constraint. 
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Figure [3] shows that the DCFs have maxima at different defays between the light curves. 
The maxima are distinct only for the optical light curves, and indicate that variations in 
B band lead variations in R band by ~0.5-l hr, while variations in V and / bands are 
simultaneous within 0.5 hr with the variability in R band. Although the peaks of the DCFs 
for the near-IR light curves are broad and their positions are uncertain within 4-5 hrs, 
the delay of variations in the near-IR region compared to those in R band is statistically 
significant. The peaks of the DCF and their delays for each light curve with respect to the 
R band light curve are listed in Table HI Such delays between the light curves are expected 
if the sizes of the emission regions at different wavelengths scale with the timescales of the 
variability found above. 



3.3. Polarization Properties 

Table 3 shows that very high degrees of polarization (up to 30%) were observed during 
the December 2006 outburst. Figure H] demonstrates that the degree of polarization increased 
when the object brightened. The latter is confirmed by the high value of the DCF peak 
between the flux in R band and the degree of polarization (Table HJ). The DCF indicates 
also that the R band and fractional polarization variations were simultaneous within 0.5 hr 
(Fig. E]). Figure H] shows a wide range, ~90°, of variability of the p olarization position 



angle. However, 9 aligns closely with the parsec jet direction, ~ —15° OJorstad et al.l 12001 



Piner et al. 1120061 ). when the degree of polarization exceeds 20%. 



Visual inspection of Figure [2] reveals a correlation between Stokes parameters - between 
Q and I and between U and I - and different behavior of the Stokes parameters from one 
epoch to another. The correlation coefficients calculated separately for each Julian date, 
given in Table [5], confirm the tight correlation between the Stokes parameters. 



4. Properties of the Source Responsible for Variability 

4.1. Method of Analysis 

Our current understanding of active galactic nuclei pictures an AGN as a complex system 
in which different components - the host galaxy, accretion disk, corona above the disk, broad 
line regions, and jet - contribute to the total flux. In some cases multicolor photometric and 
polarimetric observations of variability may make it possible to derive the characteristics of 
the variable source (spectral energy distribution, SED, and polarization parameters) without 
knowledge of the absolute contribution of the other components to the total radiation. 
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Let us suppose that the variability within some time interval is due to a single variable 
component. If the variability is caused only by its flux variation but the relative Stokes 
parameters q and u (for polarimetry) or relative SED (for photometry) remain unchanged, 
then in the space of absolute Stokes parameters {/, Q, U} or in the n-dimensional flux space 
{Fi, F n } (n is the number of spectral bands used in multicolor observations) the ob- 
servational points must lie on straight lines. The slopes of these lines are relative Stokes 
parameters of the variable component (for polarimetry) or flux ratios for different pairs of 
bands as determined by its SED (for photometry). 

With small caveats, the opposite is also true. In the case of photometry, a linear relation 
between observed fluxes in two different wavelengths during some period of flux variability 
implies that the flux ratio does not change. Such a relation for several bands indicates that 
the relative SED of the variable component remains steady and can be derived from the slopes 
of the lines. In the case of polarimetry, a linear relation between observed Stokes parameters 
suggests that the variable component has constant degree of polarization and polarization 
position angle, the values of which can be obtained from the relative Stokes parameters q 
and u that correspond to the slope of the Q vs. I and U vs. I lines, respectively. Further 



details of the method of analysis are provided in lHagen- Thorn fc Marchenko I (119991 ) . 



4.2. Polarization Parameters of the Variable Component 



Correlation plots of the absolute Stokes parameters for three consecutive epochs during 
the outburst (JD 2454077, 2454078, and 2454079) are given in Figure [5] (a,b, and c, respec- 
tively). One can see that in all cases the points lie on straight lines that leads us to conclude 
that during a given epoch the relative Stokes parameters of the variable component remained 
unchanged. The slopes of the straight lines in Figure [5] are determined by the orthogonal 
regression method. They give the relative Stokes parameters q var and u var of the variable 
component that we use to calculate its degree and position angle of polarization, p var and 
Q var (Table E]). Table [5] indicates that the degree of polarization of the variable component 
is very high and decreases from 50% to 30% as the object fades. This is accompaned by a 
systematic variation of the polarization direction. At the highest degree of polarization the 
EVPA tends to align to within 1 5° with the direction of the parsec-sca le jet inferred from 
the VLBI maps of AO 0235+164 Jjorstad et al.ll200ll : IPiner et al-fcod ). 
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4.3. Spectral Energy Distribution of the Variable Component 

We compare the fluxes in different bands relative to the R and J fluxes for the optical 
and near-IR data, respectively (see Fig. [6]), since these two bands have the largest number 
of observations. According to Figure [2J the data reveal intraday variability that can affect 
the flux-flux comparison because the measurements are not completely simultaneous. For 
example, during Julian dates with intensive monitoring, delays of ~ 10 — 15 minutes between 
observations in different bands occur. For such epochs the fluxes in R band were interpolated 
to the corresponding times of the B, V, and I bands measurements (the same for J band 
with respect to H and K bands). For epochs with a few measurements, the average at each 
band was used. 

The slopes of the lines were derived using the orthogonal regr ession method . The y 



are corrected for extinction based on the information provided by iRaiteri et al. I (120051 ). 
The logarithms of the corrected values are listed in Table HI Using these values, we have 
constructed SEDs in the optical and near-IR regions as well as the combined SED (via the 
flux-flux diagram for R and J) for the whole range from B to K bands (Fig. [7]). Figure [7] 
shows that the spectra can be represented by a power law, F v ~ v~ a , where a is the spectral 
index. In the near-IR region a = 1.10 ± 0.03, while in the optical region a = 1.39 ± 0.12. 
Although there appears to be a slight steepening of the spectrum in the optical region, this 
falls within the 2a uncertainty, so that the combined SED can be fit by a single power law 
with a = 1.04 ±0.06. This indicates that the same component with the relative SED shown 
in Figure [7] is responsible for the source variability in the whole range from B to K bands. 



5. Discussion 



The results of the polarimetric and photometric observations shows that the dramatic 
variability seen in AO 0235+164 in 2006 December might be explained by a single variable 
component that possesses the following properties: (i) steady spectral shape; (ii) increase 
of size of the emitting region with wavelength; (iii) very high degree of polarization, up to 
50%; (iv) strong correlation between degree of polarization and flux; and (v) tendency of 
the magnetic field at the highest polarization state to be perpendicular to the inner radio 
jet direction. These properties are characteristics of the model that describes an outburst in 
a bla zar as a transverse shock propagating down a relativistic jet with tur bulent magnetic 
field (IMarscher fc Gear Ill985l ; iHughes. Aller. fc Aller Ill985l ; iMarscher Ill996l ). The timescale 
of variability relates to the thickness of the shock front, which is determined by the lifetime 
of relativistic electrons accelerated at the front. This thickness decreases with frequency 
owing to synchrotron radiation losses. The shock orders the turbulent magnetic field along 
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the front, which is perpendicular to the jet direction. 

We consider a model of a square-wave transverse shock moving through a turbulent 
plasma with a constant Lorentz factor, T. The size and shape of the emission region at a 
given frequency are steady. In the observer's frame the flux of the shocked region is boosted 
as follows: 

F = F z/- Q 5' (2+Q) 5 (3+a) , (1) 

where 5' is the Doppler factor of the shocked plasma in the rest frame of the shock front 
and S is the Doppler factor of the shock front in the observer's frame. Determination of 5' 
requires knowledge of the speed and equation of state of the plasma behind the shock front 



(ICawthorne Ill99ll ). which is difficult to specify due to a lack of the necessary information. 
Instead, we assume that the velocity of the shocked plasma in the frame of the shock front 
<C c and 5' ~1. We can then approximate that 5 = [T(l — (3 cos$)] _1 , where (3 = \/I — T~ 2 
is the speed of the shock and $ is the angle between the line of sight and the velocity of the 
centroid of the emission region. 

The degree of polarization depends on the viewing angle of the shock, and ratio 
of the density in the sh ocked region with respect to unshocked region, r\ = ^shock/^unshock 



flHuehes & Miller Ill99lh : 



a + 1 (1 - r)~ 2 ) sin 2 * 
a + 5/3 2- (1 -rr 2 )sin 2 ^' 



The viewing angle of the shock in the observer's frame is subjected to relativistic aberration 
and determined by the bulk Lorentz factor and viewing angle of the jet: 



^ = tan _1 [sin$/(r(cos$ - Vl - r~ 2 ))] (3) 

We use the derived spectral index of the variable component a = 1.04 and assume 
r = 30, consistent with the highest rel iable apparent speed seen in the jet of AO 0235+16 



( IJorstad et al.ll200ll ; iPiner et al. II2006I ). We apply equation ([T]) to the photometric data in 
the R band to estimate the Doppler factor as a function of time. The scaling factor F is 
determined under the assumption that F = F max z/* / '5i 3+a \ where F max is the maximum 
observed flux in R band and S is obtained for $ Q , which is calculated from equations (T5]) 
and (J3j) for the maximum derived degree of polarization of the variable component (Table 
|5]). In this case we adopt 7/=2.5 that is the minimum possible compression that produces 
polarization as high as 50%. This yields S = 23.3 at the maximum polarization. Figure [8] 
shows the values of the Doppler factor and viewing angle of the jet derived from the observed 
flux variability. The latter allows us to estimate the viewing angle of the shock (eq. [3]) and 
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value of the plasma compression as functions of time using the light curve of fractional 
polarization (eq. [2]). The values of $ and rj are shown in Figure [8] as well. 

Figure [8] demonstrates that the observed flux and polarization behavior can be explained 
by a small change in the viewing angle < 1° and by variations of shocked plasma compression 
while the bulk Lorentz factor remains constant. A small viewing angle (~ 2.5°) agrees with 
the highly core-dominated radio images of the blazar. A change in the viewing angle can be 
caused either by curvature of the jet by < 1° or by a curved trajectory of the shock-causing 
perturbation inside the jet for an opening angle > 1°. In the brightest flux state the viewing 
angle is the nearest to the line of sight, the shock is seen almost edge-on (\& ~ 103°), the 
plasma compression is the highest (rj ~ 1.7), and the parameters change on timescale of 
hours. In the end of the outburst rj ~ 1.1 and the parameters are steady within a week. 
Changes in the rate of injection of relativistic electrons could be related to variations in 77. 

In this scenario a shock is responsible for the injection of the relativistic electrons into 
the emitting region, leading to a steady spectral shape and a short timescale of the variability 
in the optical and IR bands. Therefore, we assume that the total flux and linear polarization 
variability is caused by changes in the viewing angle and in injection of relativistic electrons 
on timescales of hours. Radiative losses determine the size of the emitting region behind 
the shock front, which is of the order of 2-3 light hours (Table H]). A steady spectral shape 
implies a quasi-steady state of the emission, which is re-established every 2-3 hrs, between 
the rate of injection and the radiative losses. In the observer's frame the synchrotron lifetime 
of electrons emitting at frequency v is: 



t loss « 4.75 x 10 2 



1 + z 
bv GHz B% 



1/2 

days, (4) 



where Bq is the magnetic field in Gauss. Since t\ oss « t u , for 5 = 23.3 (the maximum derived 
Doppler factor) and r u = 1.7 hr in the R band (Table BJ, eq. (HI) yields an estimate of the 
magnetic field, B ^ 0.5 G, that agrees with other estimates of the typical magnetic field in 



a blazar core (e.g., iMarscher &: Gear Ill985l ). Equation (T4]) implies also that the timescale 
of variability in K band should be longer than that in R band by a factor of 1.8, consistent 
with the values of r u derived from the observations. The highest degree of polarization, 
Pvar = 50%, is close to the maximum possi ble polarization (7 5%) for a synchrotron source 



with a = 1.04 and a uniform magnetic field (IPacholczyk Ill970l ). This indicates an extremely 



well ordered magnetic field behind the shock front along the plane transverse to the jet 
direction. Therefore, we consider that the polarimetric and multicolor photometric behavior 
of AO 0235+164 during the outbust in 2006 December is consistent with the properties of 
the shock-in-jet model. 
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6. Conclusion 

In 2006 December the BL Lac object AO 0235+164 underwent a spectacular outburst 
with an increase in optical flux by a factor of 4 and with degree of polarization reaching 
30%. Analysis of the multicolor photometric data and polarization measurements in R band 
indicate that: 

(i) The flux increase has an achromatic character. 

(ii) The spectral energy distribution of the variable component is a power law with 
a=1.04±0.06. 

(iii) Variations on timescales of hours are present at all wavelengths. The minimum 
timescale of variability increases with wavelength. 

(iv) There is a delay between the variations at different wavelengths, with fluctuations 
occurring first in B band and last at the near-IR wavelengths. 

(v) The fractional polarization correlates with the flux. 

(vi) The direction of polarization at the highest polarization states tends to align with 
the inner radio jet direction. 

These properties of the variable component responsible for the outburst match the 
characteristics of a transverse shock propagating down the jet. We find that the observed 
photometric and polarimetric variability can be reproduced in detail by the shock propa- 
gating down a twisted jet such that the viewing angle changes by < 1° and the plasma 
compression decreases by a factor of ~1.5 while the bulk Lorentz factor, V = 30, remains 
constant. 

This work is supported by the Russian Fund of Basic Research, grant 05-02-17562. S. 
Jorstad and A. Marscher acknowledge support from the National Science Foundation under 
grant no. AST-0406865. 

REFERENCES 

Burbidge, G. R., Jones, T. W., & O'Dell, S. L. 1974, ApJ, 193, 43 

Carini, M. 2006, in Blazar Variability Workshop II: Entering the GLAST Era, ed. H.R. 
Miller, K. Marshall, J.R. Webb, & M.F. Aller, ASP Conference Series, 350, 55 



- 11 - 



Cawthorne, T. V. 1991, in Beams and Jets in Astrophysics, ed. P. A. Hughes, Cambridge 
Astrophysics Series, (Cambridge Univ. Press), 187 

D'Arcangelo, F. D. et al. 2007, ApJ, 659, L107 

Edelson, R. A., & Krolik, J. H. 1988, ApJ, 333, 646 

Gabuzda, D. C, Rastorgueva, E. A., Smith, P. S., & O'Sullivan, S. P. 2006, MNRAS, 369, 
1596 

Hagen-Thorn, V. A. 1972, Astron. Circ, 714, 5 (in Russian) 
Hagen-Thorn, V. A. & Marchenko, S. G. 1999, Baltic Astron., 8, 575 

Hagen-Thorn, V. A., Marchenko, S. G., Yurchenko, A. V., Yakovleva, V. A. 1994, Astron. 
Reports, 38, 635 

Hagen-Thorn, V. A., et al. 2006, Astron. Reports, 50, 458 

Hagen-Thorn, V. A., et al. 2007, Astron. Reports, in press 

Hughes, P. A., Aller, H. D., & Aller, M. F. 1985, ApJ, 298, 301 

Hughes, P. A., Aller, H. D., & Aller, M. F. 1989, ApJ, 341, 68 

Hughes, P. A., & Miller, L. 1991, in Beams and Jets in Astrophysics, ed. P. A. Hughes, 
Cambridge Astrophysics Series, (Cambridge Univ. Press), 1 

Jorstad, S. C, et al. 2001, ApJS, 134, 181 

Jorstad, S. G., et al. 2007, A J, 134, 799 

Junkkarinen, V. T., et al. 2004, ApJ, 614, 658 

Lister, M. L., & Smith, P. S. 2000, ApJ, 541, 66 

Marscher, A. P., & Gear, W. K. 1985, ApJ, 298, 114 

Marscher, A. P. 1996, in Blazar Continuum Variability, ed. H.R. Miller, J.R. Webb, & J.C. 
Noble, ASP Conference Series, 110, 248 

Mead, A. R. J., et al. 1990, A&AS, 83, 183 

Miller, H. R., & Noble, J. C. 1996, in Blazar Continuum Variability, ed. H.R. Miller, J.R. 
Webb, & J.C. Noble, ASP Conference Series, 110, 17 



Pacholczyk, A. G., 1970, Radio Astrophysics (W.H. Freeman and Company, San Francisco) 
Piner, B. G., Bhattarai, D., Edwards, P. G., & Jones, D. L. 2006, ApJ, 640, 196 
Raiteri, C. M., et al. 2001, A&A, 377, 396 
Raiteri, C. M., et al. 2005, A&A, 438, 39 
Raiteri, C. M., et al. 2006, A&A, 459, 731 

Smith, P. S. 1996, in Blazar Continuum Variability, ed. H.R. Miller, J.R. Webb, & J.C. 
Noble, ASP Conference Series, 110, 135 



This preprint was prepared with the AAS IATgX macros v5.2. 



-13- 



Table 1. Optical Data 


MJD 

(1) 


B 

(2) 


MJD V MJD R MJD 

(3) (4) (5) (6) (7) 


I 

(8) 


27.3865 
28.3267 
29.3514 


20.26±0.17 
19.71±0.17 
19.95±0.17 


27.3946 19.06±0.05 27.3784 18.43±0.03 27.3823 
28.3345 19.01±0.07 28.3489 18.44±0.04 28.3473 
29.3592 19.08±0.05 29.3749 18.45±0.04 29.3735 


17.65±0.03 
17.63±0.03 
17.68±0.03 


Note. - 


MJD is JD-2454000; the table appears in full form electronically. 








Table 2. Near-IR Data 




MJD 

(1) 


J 

(2) 


MJD H MJD K 

(3) (4) (5) (6) 




979.3865 
980.3267 
985.3514 


17.24±0.11 
17.25±0.07 
17.27±0.06 


979.5627 16.45±0.15 979.5666 15.31±0.19 
980.6206 16.32±0.09 980.6258 15.47±0.16 
985.5778 16.19±0.05 985.5828 15.70±0.16 




Note. - 


MJD is JD-2453000; the table appears in full form electronically. 








Table 3. Polarization Data in R band 




MJD 

(1) 


P,% 
(2) 


e° 

(3) 




77.3307 
77.4282 
77.4433 


32.7±1.0 149.5±0.8 
27.7±1.4 152.2±1.3 
25.9±1.6 153.8±1.8 





Note. — Columns are as follows: (1) - JD-2454000; (2) - degree of the polarization; (3) - position angle of 
the polarization; the table appears in full form electronically. 
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Table 4. Multicolor Properties of the Variable Component 



Band 




N r min ,hr 


DCF pcak 


At,hv 


ig(F v /F R ) 


(1) 


(2) 


(3) (4) 


(5) 


(6) 


(7) 



B 


14.833 


24 


(5-27t^) 


0.89±0.10 


+i.ot;- 5 


-0.211±0.012 


V 


14.736 


24 


(6.1518-.S2) 


0.95±0.11 


o-oio.5 


-0.120±0.006 


R 


14.670 


66 


i 7 o+0.22 
lm ' °-0.29 


1.0 


0.0 


0.0 


I 


14.574 


24 


(4-8i±8:S) 


0.85±0.09 


O-Oto.5 


0.142±0.007 


J 


14.387 


55 


2.911812 


0.98±0.02 


-3-5i^ 


0.256±0.007 


H 


14.262 


55 


O 1Q+055 
z ' iy -0.37 


0.98±0.03 




0.400±0.006 


K 


14.140 


54 


o ni+1-21 
O.U1_ 67 


0.97±0.03 


-3-5^ 


0.526±0.005 


P 


14.670 


54 


31 + 04 

U - J1 -0.04 


0.74±0.08 


0-011.5 





Note. — Columns are as follows: (1) - filter of photometry, P corresponds to the polarization data; (2) 
- logarithm of effective frequency of the filter; (3) - number of observations; (4) - minimum timescale of 
variability, unreliable values are enclosed in brackets; (5) - peak of discrete cross-correlation function of the 
light curve at frequency v with respect to R band light curve; in the case of polarization data the peak is 
given for the DCF between the i?-band and fractional polarization light curves; (6) - delay of the peak; 
+ means that the light curve at frequency v leads the R band light curve, while — is the opposite case; 
(7) - logarithm of the contribution of the variable component in the SED at frequency v relative to the 
contribution in R band. 



Table 5. Polarization Properties of the Variable Component 



MJD 


r QI 


rui 


Pvar,% 


P) 

^var j 


(1) 


(2) 


(3) 


(4) 


(5) 


77.447 


0.93 


-0.90 


50.0±9.4 


-30.4±5.4 


78.373 


-0.48 


-0.98 


36.7±3.4 


-49.8±2.7 


79.320 


-0.58 


-0.92 


33.1±5.2 


-53.9±4.5 



Note. — Columns are as follows: (1) - Julian date minus 2454000; (2) - correlation coefficient between 
Stokes parameters Q and I; (3) - correlation coefficient between Stokes parameters U and I; (4) - degree of 
polarization; (5) - position angle of polarization. 
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Fig. 1. 



Light curves of AO 0235+16 in the optical and IR bands. 
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Fig. 2. — Simultaneous photometric and polarimetric data for AO 0235+16 at three consec- 
utive Julian dates during the outburst. 
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Fig. 3. — Discrete cross-correlation function between light curves in different bands with 
respect to the i?-band light curve; a negative delay means that R band variations lead. 
Bottom panel: DCF between degree of polarization and total flux density in R band. 
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Fig. 4. — Left panel: Dependence between the degree of polarization and flux in R band. 
Right panel: Dependence between the degree of polarization and polarization position angle; 
the dotted line indicates direction of the parsec scale jet. 
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Fig. 5. — Comparison of absolute Stokes parameters for three consecutive Julian dates during 
the outburst (a - 2454077, b - 2454078, c - 2454079); crosses correspond to the uncertainties 
of the parameters. 
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Fig. 6. — Flux-flux diagrams for optical (a) and (b) IR wavelengths. 
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Fig. 7. — The relative SEDs at IR and optical wavelengths (dashed lines) and for the com- 
bined spectrum (solid line). The optical SED is relative to the R band flux. The IR SED 
is relative to the J band flux. The combined SED is constructed using a linear dependence 
between simultaneous observations in J and R bands. 
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Fig. 8. — The derived values of Doppler factor, 5, angle between the jet axis and line of 
sight, $, viewing angle of the shock, and compession factor of the shocked plasma, r\. 



